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Abstract 

We prove that in the limit where the (3 function is dominated by the 1-loop contribution 
("large (5q limit") diagonal Pade Approximants (PA's) of perturbative series become exactly 
renormalization scale (RS) independent. This symmetry suggest that diagonal PA's are 
resumming correctly contributions from higher order diagrams that are responsible for the 
renormalization of the coupling-constant. Non-diagonal PA's are not exactly invariant, 
but generally reduce the RS dependence as compared to partial-sums. In physical cases, 
higher-order corrections in the f3 function break the symmetry softly, introducing a small 
scale and scheme dependence. We also compare the Pade resummation with the BLM 
method. We find that in the large- Nf limit using the BLM scale is identical to resumming 
the series by a x[0/n] non-diagonal PA. 



1 Introduction 



Pade Approximants (PA's) have proven to be useful in many physics applications, 
including condensed-matter problems and quantum field theory JO. We denote Pade 
Approximants (PA's) to a perturbative QCD series, describing a generic effective- 
charge S(x) = x(l + r\x + r 2 x 2 + • • • + r n x n ) by 

x[N/M] = x lll^XlZ- ■ X[N ' M] = S + X °^ N+M+1 ^ W 

i.e. the PA's are constructed so that their Taylor expansion up to and including order 
N+M = n is identical to the original series PI. PA's may be used either to predict 
the next term in some perturbative series, called a Pade Approximant prediction 
(PAP), or to estimate the sum of the entire series, called Pade Summation. The 
reasons for the success of PA's in these different applications have not always been 
apparent. Indeed, rational functions are very flexible, and hence they are good 
candidates to approximate other unknown functions, but this does not seem as a 
sufficient explanation. In this paper we give for the first time an argument, based on 
renormalization group analysis, why PA's are especially well suited for summation 
of series describing to observables in QFT. 

Among the areas in which PA's have had remarkable successes has been per- 
turbative QCD 0|| where PA's applied to low-order perturbative series have been 
shown to 'postdict' accurately known higher-order terms, and also used to make esti- 
mates of even higher-order unknown terms that agree with independent predictions 
[0] based on the Principle of Minimal Sensitivity (PMS) || and Effective Charge 
(ECH) § techniques. 

In recent papers we used the Bjorken sum rule for deep-inelastic scattering of 
polarized electrons on polarized nucleons as a testing ground and showcase for the use 
of PA's. We showed that applying the appropriate PA to the Bjorken Sum Rule, the 
renormalization scale (RS) and scheme dependence are significantly reduced. This 
observation provided circumstantial evidence that PA are resumming correctly the 
most important part of higher-order perturbative corrections. However, in absence 
of an explicit mathematical argument, we could not show that this success was not 
coincidental. 

In this paper we provide the missing mathematical argument and prove that in 
the large (3q limit (cf. Eq. (|]) and @) when only 1-loop renormalization of the 
coupling is important (e.g. large-iV/ QCD) diagonal PA are exactly RS invariant, 
giving the same result, regardless of the RS. This results directly from the fact 
that in the large /?o limit the RS transformation of the coupling reduces to a homo- 
graphic transformation of the Pade argument. Diagonal PA's are invariant under 

*We use this notation, with one power of x out of the brackets so that S{x) would fit the usual 
notation of effective charge. Note that the diagonal case in this notation is x[N/N + 1]. 
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such transformations. Non-diagonal PA's are not totally invariant, but we show that 
they reduce the RS dependence significantly. 

In general the (3 function includes higher-order perturbative corrections which 
alter the running of the coupling- const ant with respect to the 1-loop evolution. Be- 
cause of this, PA's are not exactly invariant under the RS transformation. However, 
since in QCD with Nf = 3, 4 or 5, the 1-loop running of the coupling is dominant, 
PA's are still almost RS invariant. 

The physical significance of the (approximate) RS invariance of PA's is clear: a 
part of the contribution of the unknown high orders in a perturbative series is due 
to diagrams which renormalize the coupling-constant. The numerical importance of 
these terms is reflected in the RS dependence of the partial-sum. Therefore the fact 
that PA's significantly reduce the RS dependence implies that they correctly re-sum 
these terms. In the large (3q limit (for example, large Nf QCD) the arbitrariness in 
setting the RS can be interpreted as the freedom to replace the gluon propagator 
by one that is corrected by an arbitrary number of 1-loop insertions. These include 
both fermion bubbles, and gluon (and ghost) loops. The fact that diagonal PA's 
become exactly RS invariant in this limit, suggests that they provide an optimal 
resummation for this type of higher order corrections (without actually calculating 
them). Indeed we empirically found that PA's succeed in predicting higher order 
coefficients [l"2|, r3"| , |i"4"] in the large Nf limit. 



An alternative method to resum higher-order terms that renormalize the coupling 
constant is the BLM |J method, where one uses the fact that the leading terms in 
/3q can be identified from the leading powers in Nf. One then sets the scale of the 
leading-order term in the series so that all the known higher-order terms that are 
leading in f3 are absorbed. The hope is that also unknown higher-order terms that 
are leading in (3q will be accounted for by this choice of scale. Thus we see that at 
least in the large Nf limit, PA's and BLM are very close. Indeed we found that in 
this limit choosing the BLM scale is mathematically identical to using non-diagonal 
x[0/n] PA's. This equivalence does not hold in the general Nf case. 

A detailed comparison Q of PA's with other methods aimed at optimizing per- 
turbative expansions by setting the scale and scheme, such as ECH and PMS has 
shown that there are no algebraic relations between these two methods and PA's. 
On the other hand, good numerical agreement was found in the case of the Bjorken 
Sum Rule. In this paper we show that this agreement is general, and related to the 
approximate RS invariance of PA's. 

The paper is organized as follows: In section 2 we introduce the notation and 
discuss RS dependence in the large-/?o limit. In section 3 we prove the exact invari- 
ance of diagonal PA's under RS transformations in this limit, and demonstrate it 
by a simple example. In section 4 we discuss non-diagonal PA's and show that their 
RS dependence is reduced as compared to the corresponding partial-sums. Section 
5 is devoted to the application of our results to the physical regime of QCD with 
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Nf = 3, 4 or 5 flavors. In section 6 we present our physical interpretation of the 
results and compare the PA method to the ECH, PMS and BLM methods. Section 
7 includes our conclusions. 

2 RS dependence of observables and the large j3o 
limit 

At any order of perturbation theory there remains a residual scale dependence of 
the order of the next term in the series. Thus, as one goes to higher orders, one 
expects the scale dependence to decrease. Unfortunately, in practice observables 
can only be calculated up to some finite (and usually low) order, and thus the theo- 
retical predictions are ambiguous. In QCD observables have so far been calculated 
only up to next-to-next-to- leading order (at most). Since the coupling-constant in 
experimentally interesting energies is quite large, this theoretical difficulty becomes 
a serious problem in comparison of theory with experiment. 

The residual RS dependence at a given order is an indication of the importance 
of higher order corrections, in the sense that they should compensate for this de- 
pendence, and therefore their contribution cannot be smaller than the ambiguity 
introduced by the choice of scale. The actual situation is more complicated, since 
the ambiguity in the choice of scale is not well-defined: how far a RS can be from 
the physical scale without introducing new physics? In what renormalization scheme 
would we define it? and so on. We will not deal here with these well-known open 
questions any further, but we will make strong use of the RS dependence arguments 
mentioned above. 

Suppose we start with an effective chargef] of an observable S in some renor- 
malizable QFT, calculated up to some order n + 1 in perturbation theory in some 
renormalization scheme and scale, 

S = x (l + r\X + r 2 x 2 + r 3 x 3 + • ■ ■ + r n x n ^j (2) 

In general, both the coefficients r, and the coupling-constant | i = a s /7r depend 
on the renormalization scheme and scale. If renormalization is self-consistent then 
this dependence cancels amongst the different terms in the series, so that the total 
change in S when changing the scheme or the scale is of order x 0(x n+1 ). 

More technically, the coupling-constant x in (|2]) obeys the Renormalization 

tFor simplicity we do not deal here with the more general case of S(x) — x p (l + r\x + ■ ■ •) 
for p 7^ 1. The generalization of out results is, however, straightforward. 

'■The notation we use is suitable for QCD, and in some cases we explicitly state QCD results, 
but the conclusions apply to other QFT's as well. 
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Group (RG) equation: 



^ = P X 2 + AX 3 + P 2 X 4 + ■ ■ ■ EE /3(X) (3) 

where the first two coefficients are the same in any renormalization scheme, while 
higher coefficients, {3%, ft^, ... are renormalization scheme dependent. In fact, as 
was shown in , at any given order in perturbation theory different renormalization 
schemes can be uniquely defined by these higher coefficients of the ft function. 

In QCD the first three coefficients of the ft function are known Ul7,18,19|. The 
first two are: 

* = -\(t c --I n >) (4) 

where in SU(N C ) gauge theory C a = N c and C f = (N 2 - 1)/{2N C ). 

We are interested in studying the effect of changing the renormalization scale. 
Suppose that x above refers to the coupling-constant being renormalized at the 
physical scale (for instance, the momentum transfer Q 2 in deep inelastic scatter- 
ing), while we denote by y another legitimate coupling-constant renormalized at 
some arbitrary scale /z 2 . We denote by t the scale shift that corresponds to the 
transformation from x to y: 

t = in g) (6) 

Equation (|3|) can be integrated order by order to give the required scale-shift trans- 
formation: 

x = y + ft ty 2 + {ftlt 2 + P^y 3 + ^ 3 t 3 + + fat) y 4 + ■ ■ ■ (7) 



We turn now to the large fto limit, defined by the condition 

Po > PiX 1 (8) 

for any % > 1. This is the limit in which our main argument is given. Condition (H) 
implies that only ftfi t k terms in equation (|^) are significant. An important example 
where this condition is realized is QCD with a very large number of flavours (see 
[ |TjJ] and the Appendix). Since Po is relatively large also in other cases, including in 
QCD with a few light flavours, this approximation is worth studying. We return to 
the validity of this approximation in QCD in a later section. 

It is important to point out that the question of scheme dependence does not 
appear in the large fto limit, as all the high order coefficients of the ft function are 
not important. Thus the only remaining ambiguity is due the arbitrariness of the 
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renormalization scale. Therefore, the study of the RS dependence is much simplified 
in this limit. 



When all the higher order terms in the (3 function are neglected, the scale shift 
transformation (^) simplifies to a geometrical series: 

x = y (l + (3 ty + (3 2 t 2 y 2 + (3 3 t 3 y 3 + • • •) (9) 
which can be summed to all orders: 

x = t\t ^ 

1 - p ty 

In the following we briefly review some properties of this scale shift transformation. 

We start by recalling the fact that the transformation (|HJ) admits the following 
mathematical relation, which we refer to as additivity of scale-shifts: 

£3 = ^1+^2 (11) 

where 

x — ^ y z (12) 
x — > z 



This property implies that the set of all transformations (|Kj) with different scale- 
shifts t forms a group. It is easy to show the existence of a unit operator (t — 0), 
the existence of an inverse transformation (t < — > —t) and associativity. As for 
the question of whether the set is a closed one: from the mathematical point of 
view - it's clearly not closed, due the the existence of the Landau pole, where the 
denominator of ( |T0D is zero: t Landau = From the physical point of view we can 
assume that the Landau pole is not reached and thus ignore this problem, concluding 
that RS shifts are indeed a group. 

As we will soon be interested in applying transformations like (|10| ) to partial sums 
and to PA's, it is important to note here that the transformation ( |I0D is a rational 
polynomial rather than a polynomial of finite order, and as such, it transforms a 
polynomial of finite order into a rational polynomial, but transform an a rational 
polynomial into another rational polynomial. 

This is a crucial observation which is the basis of the exact RS invariance of 
diagonal PA's in the large (3o limit, and the reduced RS dependence of non-diagonal 
PA's in this limit. It is also the reason why partial sums at a given (finite) order 
can never be RS invariant. 
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3 RS invariance of diagonal PA's in the large (3q 
limit 



The statement we prove is that in the large (3$ limit, starting with a partial sum 
S(x) of a given order n + 1, as in eq. (|2]) f\, the diagonal x[N/N + 1] PA of S(x), 
such that n = 2N + 1, does not depend on the RS in which the given partial sum 
was calculated. 

The following diagram will be helpful in illustrating our discussion: 



partial-sum: 


S(x) i 


— ► Si{y,t) 


— ► S 2 (x,t) 




I 


I 




PA: 


P(x) 


Pi(y,t) 


P 2 (x) 




1 


I 


1 


Taylor of PA: 


T(x) 


Ti(y,t) 





where ^(x) is the (n + l)-th order partial sum as in eq. (0), P(x) is the x[N/N + 1] 
PA constructed from S(x), and T(x) is the infinite-order Taylor series of P(x). The 
horizontal arrows represent the application of scale-shift transformation according 
to equation (|ToD , once from x to y (with the scale-shift t) , and then back, from 
y to x (— t). The first scale-shift in the first line is intentionally represented by a 
different arrow (i — >) than the other scale shifts ( — >), meaning that the series in y is 
truncated at the (n+l)-th order after applying the scale-shift transformation Q This 
truncation does not apply to the other scale-shift transformations in the diagram, 
where an exact transformation according to (|H]) is meant. In the following we 
discuss the elements of the above diagram in detail. 

The starting point in a resummation process is an (n + l)-th order partial- 
sum which generically can be either S(x) or S\(y,t). The two partial sums give 
different numerical results for the observable, but a priori one is just as good as 
the other. The numerical difference between S(x) and S\(y,t) can be calculated 
by applying the full inverse scale-shift transformation (—t) to S\(y,t) which yields 
5*2(2;, t). The latter is numerically the same as S\(y,t) but differs from the original 
S(x) by corrections of order xO(x n+1 ). 

§ We take n to be an odd number, in order to construct a diagonal PA. 

^We remind the reader that in x[N/N + 1] PA we refer to x times a rational polynomial with 
numerator of order N and a denominator of order iV + 1, as defined in the Introduction. Therefore 
we call it a diagonal PA. 

* In order to calculate Si(y), one substitutes a; as a function of y in S(x) according to (pp|), and 
Taylor expands the result to its (n + l)-th order, neglecting all the higher order terms. 
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Having described the RS dependence of the partial-sums we now consider PA's. 
We will now show that diagonal PA's are RS invariant, i.e. i^O^) does not depend 
on t. We will explicitly show that the two PA's - P(x) that is based on S(x), and 
P2{x) which is calculated by an inverse scale-shift transformation of Pi(y,t), i.e. 
based on Si(y,t), are exactly equal! This invariance results from the fact that in 
the large-/9o limit, RS transformations of the coupling amounts to a homographic 
transformation of the Pade argument. Diagonal PA's are known to be invariant 
under such transformations (see Ref. |16| and references therein). The proof is as 
follows: 



a) P(x) is calculated as the x[N/N + 1] PA of S(x): 

l + a 1 x + ... + a N x N 

P(X) = X- : : Tf— (13) 

v ' 1 + b x x + ... + b N+1 x N+l 
Similarly, P\(y,t) is calculated as the y[N/N + 1] PA of Si(y,t). 

b) The scale-shift transformation ( |i0|) is applied to P(x), to give a function which 
we denote P*(y,t): 

P;iv , t) = ( » ) 1+a ;(^f) + - +a "^)l (u) 



-Poty) 



We shall see that P*(y,t) is equal to P±(y,t). 



c) By multiplying both the numerator and the denominator by (1 — fioty) N+1 we 

see that P*{y, t) is actually a a rational polynomial of the type y[N/N + 1]. 

d) T*(y,t) is the Taylor expansion of P*(y,t). The n first coefficients of T*(y,t) 
are necessarily the same as those in S\(y,t). This is because T*(y,t) can also 
be viewed as the scale-shifted version of T(x), which, by the way the P(x) PA 
was constructed, has the same n first coefficients as S(x). On the other hand, 
the scale-shift transformation is such that the n first coefficients of T*(y, t) (or 
Si(y,t)) depend only on the first n coefficients in T(x) (or S(x)). 

e) Since an x[N, N + l] PA is uniquely determined by the first n terms of a power 

series (n = 2N + 1), it follows that P*(y,t) equals to Pi(y,t). 

f) If we now apply an inverse scale-shift (— t) on P 1 (y,t), we get Pi{x) which is 
equal to P(x) (and therefore does not depend on t). This is due the existence 
of an exact inverse scale-shift transformation, or putting it differently, due to 
additivity of scale-shifts, discussed in the previous section. [] 



TThis additivity of scales breaks down when higher powers oft are present in the transformation, 
due to non-negligible higher-order corrections to the [3 function. 
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This completes the formal proof of the invariance of PA's in the large (3q limit, as 
stated at the beginning of the section. We now present a simple example to illustrate 
how this symmetry works in practice. We start with a fourth order partial-sum 
(n = 3) ' ' 

S(x) = x (l + ri x + r 2 x 2 + r 3 x 3 ) (15) 
and calculate P(x) as a x[l/2] PA: 

P(x) = x (r,-r,') + (2nr 8 -r,'-r t )» 

(— r 2 2 + ri r 3 ) x z + (— r 3 + r 1 r 2 J x + r 2 - ?v 

Expanding P(x) back in a Taylor series we get T(x): 



r 2 3 - 2 r 2 ri r 3 + r 3 2 

~2~ 



rr-i/ \ . 2 i 3 i 4 i Z 1 3 3 5 /0 / t> \ /-irrN 

T(x) = x + rix +r 2 x +r 3 x H x + 0{ x ) (17) 

r 2 - tV 



We identify the r 4 PAP as the coefficient of x in (|i~7]). 

We calculate the RS dependence of the partial sum S'(x) as follows. First, we 
obtain Si(y, t) by substituting x = 1 _^ oj/ in 'S'(x), expanding the result into a power 
series and truncating beyond the y 4 term: 

S 1 (y,t) = y+(0 o t + r 1 ) y 2 + (f3 2 1 2 + 2 n (3 t + r 2 ) y 3 + (18) 

+ (Po 3 t 3 + r 3 + 3r 1 p 2 t 2 + 3r 2 p t) y 4 



We now transform back to x by substituting y = 1+ ^ QX in Si(y,t), and expanding 
to all orders in x, to get S%(x, t). The resulting formula for S^x, t) is 

^(x, t) = x + ri x 2 + r 2 x 3 + r 3 x 4 + (19) 

(-4 /5 3 t 3 n - /3 4 t 4 - 6 r 2 /3 2 t 2 - 4 r 3 fat) x 5 + 0( x 6 ) 

The RS dependence of the partial sum is the difference AS(x,t) between ^(x, t) 
and S(x): 

AS(x,t) = S 2 (x,t) - S(x) (20) 
= (-4 A, V n - p H 4 - 6 r 2 /5 2 t 2 - 4 r 3 /5 t) x 5 + 0( x 6 ) 

We note that the t dependence of S 2 (x,t) and thus also AS(x,t) can be quite 
large for certain values of t (at least for large t it is clear that the contribution due 
to /3 4 t 4 cannot be canceled by other terms). On the other hand, the x[l/2] PA is 
exactly RS invariant: if we start with the fourth-order partial sum Si(y,t) of Eq. 
(|i~8D, construct a J/ [1/2] PA P\(y,t) and transform back to x, we get exactly the 
P(x) of Eq. (fHf ), as proven. 

We stress that also the PAP's are RS independent. For instance, if we calculate 
the PA-improved 5-th order partial sum in y, and transform it back to the coupling x, 
we will get exactly the same prediction for r 4 , as in (|T?D, without any t-dependence. 
This is true not only for the PAP of r 4 but for any higher-order diagonal PAP. 
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4 The reduced RS dependence of non-diagonal 
PA's in the large /3q limit 



In the previous section we have proved that diagonal x[N/N + 1] PA's, are exactly 
RS invariant in the (3 limit. In this section we will see that non-diagonal x[N/M] 
PA's are not exactly RS invariant even in this limit. However, on the global level (i.e. 
for large variations of the scale t), their RS dependence is much reduced compared 
to partial sums. 

The RS dependence of non-diagonal PA's is illustrated in the following diagram 
(cf. analogous diagram in the previous section): 

partial-sum: S(x) ^ Si(j/,*) ^ S 2 (x,t) 

i i 
PA: P(x) -U Pf(y,t) ^ P^t) 

1 1 
Taylor of PA: T(x) Ti(y,t) T 2 (x,t) 



To illustrate this, we turn now to the simplest example, where we start with a 
third order series (n = 2): 

S(x) = x ( 1 + T\ x + r 2 x 2 ^ (21) 
Calculating its a; [1/1] PA, we get 

P(x) =x n + " r2) X (22) 
n - r 2 x 

Expanding back in a Taylor series we obtain: 

T(x) = x + nx 2 + r 2 x 3 + ^x 4 + 0(x 5 ). (23) 
Applying the scale shift transformation to S(x) yields: 

.2 \ 



/ ny r 2 y 
V 1 + ~, 7TT- + 



l-fcty (l-(3 t y y 



(24) 



V 1 - tp y J 1-Poty 
which can be Taylor expanded in y, and truncated at fourth order, to give: 

S 1 (y,t) = y+ {n + f3 t) y 2 + (2 n fat + r 2 + (3 H 2 ) y 3 (25) 

9 



Constructing a y[l/l] PA, we get: 



A(y,t) = y , (rx + M + ^-r^v 

-2 n - r 2 - p t 2 )y + n + p t 



On the other hand, if we transform P(x) to y, using the exact scale-shift transfor- 
mation (0), we get: 

P*(v t) = v n + (A^IlM^k (27) 

1 ^ ) y n- (r 2 + 2r 1 /3 t)y + (r 2 + r^/W 1 J 

We see that P±(y,t) is a diagonal rational polynomial, of type y[l/2], rather than 
an off-diagonal rational polynomial of the type 2/ [1/1] , like Pi(y,t) of Eq. (EBI). 



Therefore P*(y,t) ^ Pi(y,t), and the scale invariance property does not hold here. 

In order to measure the RS dependence of the non-diagonal PA, we define 
AT(x, t), in analogy with AS(x, t) which was introduced in Eq. fl20"|) to measure the 
RS dependence of partial-sums. 

AT(x,t) = T 2 (x,t) - T(x), (28) 

where as T(x) is the Taylor expansion of P(x) and T 2 (x,t) is obtained by applying 
an inverse transformation (—t) to the PA-improved partial sum T\(y,t), the Taylor 
expansion of Pi(y,t) (see diagram). 

Going back to our example, we calculate T 2 (x,t), 



T 2 {x,t) = x 



Pot + n 



(29) 

As expected, up to r 2 we obtain the same coefficients we started with. The predicted 
r 3 turns out to be scale dependent (and therefore different from the one obtained 
from P(x)). Thus for the RS dependence of the PA function AT(x,t) we obtain: 

AT( = f -Mii±Ap^±im. _ d) x * + (A (so) 

\ Pot + r 1 rij 

This has to be compared with the RS dependence of the partial sum. Taking the 
inverse (— t) transformation of Si(y,t) we get S 2 (x,t): 

S 2 (x, t) = x + n x 2 + r 2 x 3 + (-Po 3 t 3 - 3 n P 2 t 2 - 3 r 2 p t) x 4 + 0(x 5 ) (31) 
and therefore: 

AS(x, t) = (-P 3 t 3 - 3 n Po 2 t 2 - 3 r 2 p t) x 4 + 0{ x 5 ) (32) 
Comparing AT(x,t) with AS(x,t) we make the following observations: 
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a) The asymptotic behavior of AT(x, t) at large scale variations (large t), is much 

milder than that of AS(x, t) 

AT(x,t) ~ {r\ - r 2 )(3 tx 4 + 0(x 5 ) (33) 

vs. 

AS(x,t) /3 3 tV + C(x 5 ) (34) 

b) On the other hand, there is a pole in AT(x,t). At certain scales it results in 
very large deviations of the PA from the value typical at most other scales. 
We thus have to be careful not to use the a; [1/1] PA at these scales (for any 
specific case one can plot the a; [1/1] PA as a function of the RS, and identify 
the scales that are badly influenced by the pole). 



We now turn to a second example, where we construct a a; [0/2] PA from the 
same partial sum (pip. In this case we do not give all the details, but only state the 
final result for AT(x,t): 

AT(x,t) = [-A)t(r? - r 2 ) + 2r 2 r l - rf\ x 4 + 0(x 5 ) (35) 

We notice that for the a; [0/2] PA (as for the a; [1/1] PA) the asymptotic behavior 
is mild, but (in contrast to the a; [1/1] PA) no poles appear in AT(x,t). Thus we 
conclude that the a; [0/2] PA is much less scale dependent than either the partial 
sum or the a; [1/1] PA. 

An interesting result is that the asymptotic behavior of AT(x, t) for large enough 
t in the case of the a; [0/2] PA is the same as that of the a; [1/1] PA. This can 
be confirmed by comparing the large t behavior of Eq. ( p0|) and Eq. (|35|). The 



mathematical reason for this similarity is that in both cases, the non-diagonal PA's 
transform under the scale-shift transformation into diagonal rational polynomials of 
order 2 (x[l/2] PA). Since we lack one parameter (r^) in order to write a "correct" 
a; [1/2] PA to describe the observable, we are left with some ambiguity, having a full 
set of functions of the type a;[l/2] with coefficients that depend on one free parameter 
(r^) rather than one specific a; [1/2] PA. The remaining ambiguity is reflected both 
in the freedom to choose among the a: [1/1] and a; [0/2] PA's, and in the freedom to 
set the RS. The interpretation of this result is that the choice between x[l/l] and 
a; [0/2] PA is equivalent to the choice between different RS's. We will later see how 
this generalizes to higher-orders. 

Note that throughout the analysis we did not make any specific assumptions 
about the perturbative coefficients r, of the observable under consideration. In |4[] 
we examined the [1/1] and [0/2] PA's considered here, as well as the partial sum, 
for the physical example of the Bjorken Sum Rule with Nf = 3. We found that 
indeed the RS dependence of the a; [0/2] PA is much reduced as compared to the 
partial sum. We also found that the a; [1/1] PA has a large RS dependence in specific 
RS's. We now see, in retrospective from the large @ limit, that these features are 
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general. In the next section we will show why the large (3q assumption is a good 
approximation for QCD with 3 flavours. But before doing so, we want to see how 
the present conclusions for the n = 2 case can be generalized to higher n. 

For n = 3, we start with a partial sum as in Eq. (]T5f) . We can construct two 
non-diagonal PA's: #[2/1] and #[0/3], in addition to the diagonal #[1/2] PA studied 
in detail in the previous section. We found that the partial sum has a considerable 
RS dependence, as described by (|20|), while the #[1/2] PA is exactly RS invariant. 
Here we look on AT(x, t) for the #[2/1] and a; [0/3] PA's. 

The asymptotic behavior of AT(x,t) at large scale shifts is the same for #[2/1] 
and x[0/3] PA's: 

AT(x,t) ~ {r\-r 2 )f3lt 2 x b + C(x 6 ) (36) 
to be compared with the asymptotic behavior of AS(x,t) from equation (f20|): 

AS(x,t) /3 4 t 4 x 5 + 0{x 6 ) (37) 

AT(x,t) for the #[2/1] PA has poles (the coefficient of x 5 has two poles in t), 
whereas AT(x, t) for the a; [0/3] PA does not. This makes the use of the x[0/3] PA 
much safer than the use of the #[2/1] PA. However, here it is preferable to use the 
diagonal #[1/2] PA which is exactly invariant and therefore the other two will not 
be relevant. 

For n = 4 we have four non-diagonal PA's: #[3/1], x[2/2], #[1/3] and #[0/4]. 
The asymptotic behavior of AT(x,t) for #[3/1] and x[0/4] PA's is the same, and 
the leading term is proportional to /3gt 3 x 6 ; the asymptotic behavior of AT(x,t) for 
#[2/2] and #[1/3] PA's is again the same, and the leading term is proportional to 
flotx 6 . These are to be compared to the asymptotic behavior of the RS dependence 
of the partial sum (AS(x,t)) in which the leading term is /3^t 5 x 6 . Out of the four 
PA's, only a; [0/4] does not have poles in AT(x,t). 

Generally, starting with a partial sum of order n + 1 and constructing an x[N/M] 
PA (JV + M = n), we find the asymptotic behavior of 

AT{x,t) ~ $t d x n+2 (38) 

where d is given by: 

d = \N + 1- M\ = \2N + l-n\ (39) 

while the asymptotic behavior of AS(x,t) is 

AS{x,t) ~ p™ + H n+1 x n+2 (40) 

d actually measures the dimensionality of the ambiguity space (= number of un- 
known parameters) that exists in writing a diagonal x[N/N + 1] PA (for iV > M) 
or x[M — 1/M] PA (for N < M — 1) describing the series. It's clear that d < n + 1 
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and therefore the global RS dependence of the x[N/M] PA is milder than that of 
the partial sum. On the other hand, already at the leading term in AT(x, t) there 
will be N poles which will cause specific RS's to exhibit sharp scale dependence, 
something that we want to avoid. Only the a;[0/n] PA will not have any poles in 
AT(x,t). 

We conclude that in general non-diagonal PA's also have a reduced RS depen- 
dence as compared to the partial-sum. In Q we recommended considering the RS 
dependence in choosing the appropriate PA. The above discussion shows that in 
cases where the large (3q approximation is valid, for series of even orders in the 
coupling-constant (odd n = 2N + 1), the diagonal x[N/N + 1] PA's should be pre- 
ferred, while for odd orders in the coupling (even n = 2N), the next to diagonal 
PA's (x[N/N\ and x[N — 1/N + 1]) are best from the point of view of the global 
RS dependence, while the a;[0/n] PA is the only one in which poles in AT(x,t) are 
guaranteed not to appear, and thus likely to exhibit the least RS dependence for 
small scale variations. 

5 Applicability of the results for QCD with Nf = 

3, 4 or 5 

In this section we leave the large (3q limit and consider the physical example of QCD 
with 3, 4 or 5 quark flavors. 

As mentioned in the Introduction, we already have strong evidence @,||,|J that 
PA's are useful in this case. Moreover we showed in Q that x[0/2] PA significantly 
reduces the RS dependence of the Bjorken Sum Rule in any renormalization scheme. 
Therefore we already know that the results we presented here for the large (3q limit 
are approximately true also in real- world QCD, at least for the example of the 
Bjorken sum rule. In this section we will show that the approximate invariance of 
PA's under RS transformations holds for a generic QCD observable in a wide range 
of renormalization schemes. Actually, we believe that this approximate invariance 
holds in many QFT examples, since the only requirement is that the running of the 
relevant coupling will be dominated by the 1-loop contribution. 

In general, the scale-shift transformation (|7]) cannot be written in the form (|T0D, 
since higher order terms in the (3 function do not vanish. In fact, the problem 
becomes even more complicated because the generic scale-shift transformation is 
not well defined, since it can only be written as a divergent asymptotic series, and 
might not even be Borel-sumable, due to IR renormalons. In addition to the RS 
dependence, there is also renormalization scheme dependence, which can be viewed 
as arbitrariness in setting the higher order coefficients of the (3 function. This makes 
it clear that the results of RS invariance in the large (3q limit cannot be formally 
extended to the general case. 
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On the other hand, we can check to what extent the scale-shift transformation 
which is exact in the large fa limit can serve as an approximation to the actual 
scale-shift in the general case (0). Let's consider the first two terms in (0) which 
deviate from fllCf) - these are the third and the fourth orders in y: 



f3 2 t 2 — (3 2 t 2 + fat 

ffit 3 — (31$ + \fafat 2 + fat 

In order to measure the effect of these higher-order terms we define: 



(41) 



™ = ®W* (42) 

and 

CaW ^3 ( 43 ) 

We see that for any (non-zero) fa and large enough scale shift t, the dominance of 
the leading /3gt & terms is recovered. 

Since we are interested specifically in QCD with Nf = 3, 4 and 5, we can 
numerically estimate the dominance of the leading /3^t k terms. For instance, for 
N f = 3 we have in the MS renormalization scheme 000: fa = -2.25, fa = -4 
and fa = —10.06. Using these values we calculate Ci(t) and C<z(t) and look at the 
actual transformation of the coupling as a function of t. 

In Figure [I] we plot y(t) for x = a s /n = 0.07, corresponding to Q 2 = 20 GeV 2 
in MS with fi 2 = Q 2 . We calculate y(t) in three different ways: 

a) The exact transformation in the large fa limit: 

»<«> = T+7wJ < 44 > 

b) The first four terms in the leading fa transformation. 

y (t) = x - fatx 2 + (3 2 t 2 x 3 - f3lt z x A (45) 
(the difference between (|45|) and (p|) at large t is due to higher orders). 

c) The first four terms in the QCD RS transformation, with Nf = 3 (fa taken in 

MS). 

y (t) = x - fatx 2 + d(-t)p 2 t 2 x 3 - C 2 (-t)plt 3 x 4 (46) 

Figure |2| presents the corresponding relative deviations of the the leading fa 
transformations ( 44]) and (^5| ) from the QCD transformation fl4"6f) . 



From figures |l| and ^| we see that the leading fa approximation is quite accurate 
for QCD with 3 flavours in a very large range of scales. We also find that already 
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at \t\ ^ 4 the error due to the neglected higher order terms becomes larger than the 
one due to neglecting of (3\ and /3 2 - 

The numerical results for Nf = 4,5, as well as for other renormalization schemes 
(i.e. other values of p 2 ) are almost identical to those presented in figures [l] and |[ 

The conclusion from this analysis is that the scale- shift transformation in QCD 
with Nf = 3, 4 or 5 can be well approximated by the large (3$ scale-shift trans- 
formation of ([44]) . Therefore, to a good approximation, the conclusions we drew 
concerning the invariance of PA's in the large /3q limit apply in the realistic case 
as well. As mentioned, we already found through an explicit calculation [|J, that 
a; [0/2] PA reduces the RS dependence of the Bjorken Sum Rule in wide range of 
renormalization schemes. 

6 The physical interpretation and comparison with 
ECH, PMS and BLM 

In the previous sections we showed that diagonal PA's become RS invariant in the 
limit of large flo, while non-diagonal PA's generally reduce the RS dependence of 
the observable. Since the leading /3q approximation is good in many physical cases, 
including QCD, PA's have a much reduced RS dependence there as well. In this 
section we discuss the physical interpretation of this result and compare the PA 
method to the ECH, PMS and BLM methods. 

6.1 Physical interpretation 

In general, a significant part of the contribution of unknown higher orders in a 
perturbative series is due to diagrams that renormalize the coupling-constant. The 
numerical importance of these terms is reflected in the RS dependence of the partial- 
sum since it is the higher order terms that compensate for this unphysical depen- 
dence. 

In order to understand the meaning of the (approximate) RS independence of 
PA's, we first analyze it in the large /3q limit, where the running of the coupling- 
constant is completely determined by a 1-loop renormalization. 

We found that in this limit diagonal PA's become exactly RS invariant. This 
strongly implies that the diagonal PA sums-up the higher order contributions that 
compensate for the RS dependence of the corresponding partial-sum. 

On the other hand, even in the large /3q limit, we do not expect the diagonal PA 
to sum the higher order terms exactly. This is because rational functions cannot 
represent factorial behavior of the coefficients as expected at large orders due to 
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renormalons. 



All order perturbative calculations for some observables in QED and QCD in 
the large Nf limit have been performed (see for instance [12] and |Tj|). The con- 



ditions for the dominance of the 1-loop running of the coupling-constant (J8|) hold 
in this limit (see the Appendix), and therefore our conclusions concerning the RS 
independence of diagonal PA's hold as well. Thus this limit provides a good test- 
ing ground for our method. In the large Nf limit only fermion-bubble renormalon 
chains contribute at high orders, and therefore the results are extremely simple in 
the Borel plane. The only remaining ambiguity is related to the integration through 
poles located on the positive real axis in the Borel plane, i.e. the IR renormalons in 
QCD (and UV renormalons in QED). After performing the integration and going 
back to the a s plane, we do not get a rational polynomial in a s . Therefore PA's 
in the a s plane cannot be expected to give the "exact" result, but only to provide 
good approximations, and "converge" with increasing orderfj]. 

Since the RS independence of PA's in the large (3q limit is exact, we suspect that 
diagonal PA's, although do not sum the full series, do provide some kind of optimal 
resummation of higher order corrections. We expect, for instance, that in the large 
fi Q limit diagonal PA's are more accurate than any other non-diagonal PA and at 
least as accurate as any scale-setting method. We show later that in the large Nf 
limit, choosing the BLM scale is equivalent to using the x[0/n] non-diagonal PA. 

Given the RS independence of PA's, it is interesting to examine numerically the 
convergence of increasing order PA's and the precision of the Pade Approximant 
predictions (PAP's) as compared to exact calculations in the large- Nf limit. An 
empirical study of this kind has been done in and ||. It was found that increasing 
order PA's do not converge but oscillate around the Cauchy Principal- Value of the 
inverse Borel integral. These oscillations are due to IR renormalons and in the 
absence of such, increasing order PA's converge to the Borel sum of the asymptotic 
series. It was also found that increasing order PAP's become very close to the exact 
perturbative coefficients, and that the errors decrease exponentially with order. The 
errors can be approximated by a simple function, and thus a further increase in the 
PAP precision is possible (this issue has been discussed in H and in |T5|). 

We stress that the exact RS independence of diagonal PA's (and the reduced 
dependence of non-diagonal PA's) holds whenever the 1-loop running of the coupling 
is dominant (cf. Eq. (||)). Large (3$ does not necessarily imply large Nf, as there 
may be many other cases in QFT where the condition holds. As we saw in the 
previous section, the large /3 approximation is quite good in QCD with only 3 to 5 
flavours. However, in contrast with QED, the 1-loop running of the coupling in QCD 
is dominated by gluon (and ghost) loops and not by fermion loops. The BLM M 
method resums 1-loop insertions that renormalize the coupling-constant. Therefore 

*We use quotation marks to indicate that the terms are not well defined mathematically (due 
to the existence of singularities on the integration path). 
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it is natural to compare the PA and BLM methods, as we do in Sec. [6.3]. 

We see that the basic physical idea behind PA's, and the scale and scheme setting 
methods is the same: part of the contribution of unknown higher order corrections 
is related to the running of the coupling constant. Thus it may be possible improve 
the perturbative result through resumming part of these unknown terms by relying 
on the characteristics of the renormalization group. In the following subsections we 
compare the PA method to ECH, PMS and BLM. 

6.2 Comparison of PA's with ECH and PMS 

The method of Effective Charges (ECH) || is based on choosing the renormalization 
scale and scheme, such that the series reduces to a leading order term (all the other 
known coefficients in this scheme are exactly zero). There is a unique scheme that 
fulfills this criterion. 

In H] we found good numerical agreement between x[0/2] PA of the Bjorken 
sum-rule (which turned out to be almost RS invariant), and the ECH (and PMS) 
scheme-setting methods. As shown in [|J, there seems to be no algebraic relation 
between these two methods and the PA. Still, our results in this paper imply that 
the numerical agreement is general. The reason is quite simple, and it is directly 
related to the approximate RS independence of the PA's. Suppose a PA was exactly 
scale and scheme invariant. Then we would get exactly the same numerical result 
in any scheme and scale, and in the ECH in particular. However, in the ECH 
scheme, the series reduces to a leading order term (all the higher order coefficients 
are exactly zero) and therefore the PA is identical to the "partial sum" . Therefore 
the ECH result is exactly equal to PA in any scale and scheme. Of course, this exact 
agreement breaks down as effects of RS dependence are turned on; these include the 
possible choice of a non-diagonal PA, higher order corrections to the j3 function, and 
scheme dependence. Still, we expect that at scales and schemes not too far from 
ECH, the PA result will be much closer to the ECH result than to the corresponding 
partial-sums [| 

Another approach to set the scale and scheme is the Principle of Minimal Sen- 
sitivity ||. In this method, one chooses the scheme in which the renormalization 
scheme and scale dependence vanishes exactly. PMS is close to ECH both in its 
nature, and in it's numerical predictions (see J!]] and pj). Knowing now that ECH 
and PA's methods generally agree, we also expect PMS to be close to PA's. We note 
that reducing the RS dependence is the common basis for both PA's and PMS. The 
difference is that PA's can be applied at any scale and scheme and they reduce the 
RS dependence globally (even for large scale variations), while PMS chooses a scale 
and scheme such that the local scale dependence vanishes. 

§ This is true as long as pole effects of non-diagonal PA's are avoided. 
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6.3 Comparison between PA's and BLM 



The BLM method || is based on the observation that in the large- Nf limit of QCD 
all the higher-order corrections are due to fermion loops, which are also responsible 
for the running of the coupling. Therefore they can be absorbed by changing the 
RS at which the coupling is defined. It is reasonable to expect that absorbing these 
corrections in the scale will improve the perturbative result. 

In QCD with only a few flavours, one can use the fact that (3q is linear in Nf to 
single out the 1-loop corrections to the coupling by identifying the terms that are 
leading in Nf. Thus also in QCD 1-loop corrections to the coupling can be absorbed 
by setting the scale of the leading term. Technically this is done by setting the scale 
so that at higher orders all the leading terms in Nf will cancel exactly. 

BLM was generalized to account for non-leading corrections to the renormal- 



ization of the coupling in several ways (see, for instance, [|T(J], and fl!!] ). Since 
the basic intuition of BLM and its generalizations relies on the large Nf limit, we 
present a detailed comparison with PA's for this case only. 

We start with an effective charge of a generic observable in the large Nf limit 
(see the Appendix). The BLM scale-setting procedure is based on eliminating the 
Nf dependence of the coefficients r*. In the large Nf limit this results in complete 
elimination of the r^'s, since in this limit r, oc Nf (there are no sub- leading terms), 
leading to the result: 

Sblm = x{t BLM ) (47) 
or, equivalently, using the large Nf limit notation in the Appendix: 

SSL7T = <t BLM ) (48) 

where z = xNf and 

tBLM =h + t 2 Z + t 3 Z 2 + ■■■ + tnZ™- 1 (49) 

where t\ is proportional to ri, t 2 = t 2 (ri, ^2), h = t 3 (ri,r 2 ,r 3 ), etc. 

The leading-order BLM scale t\ is chosen such that c\ (or n = C\Nf) is elim- 
inated. Using this scale results || in a summation of the leading diagrams which 
correct the gluon propagator at higher orders. For instance, a contribution like c\ z 3 
(or r\x 3 ) is accounted for (cf. ( |A.5 )), although terms of order z 3 were not initially 
included. 

In order to eliminate also the next coefficient c 2 , one has to alter the scale-shift 
by adding a term t 2 z that is proportional to the coupling^]. Similarly, one adds t 3 z 2 
to eliminate c 3 , and so on. In such a way all the known terms can be absorbed 



^ We will later discuss the effects of the non- leading corrections to tBLM- 
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into the definition of the coupling- const ant, hopefully summing correctly the bulk 
of higher-order unknown contributions. 



Several different proposals were made for generalizing BLM beyond the leading 
scale t BLM = ti, but in the large Nf limit they all agree: both the single-scale BLM 



method of Ref. || and the multi-scale method of Ref. |T0| reduce then to (f48|) and 

(EDI 

Suppose we want to calculate the effective charge in the large (3q limit by the 
BLM prescription, according to (|^) and (f^9"D, assuming the coupling-constant at 



the physical scale (i.e. z) is known. The effective charge we are calculating is just 
the coupling- const ant at the scale tBLM, given by the inverse of relation (|A.6|) : 

s - w = rr^i < 50 > 

where from now on (3 stands for to (3^^°° = |. We substitute ([49]) into (|50| ) to 
get: 

= 1 + (3 t lZ + fottZ* + ■■■ + (3 tnZ n (51) 

By construction, if S B [~^°°(z) * s expanded in powers of z up to order n + 1, 
one would get the original series S(z). We note that the r.h.s. of (pTj) is a z[0/n] 
rational polynomial. Now, since there is a unique z[0/n] PA which has an (n + l)-th 
order Taylor expansion equal to S(z), we conclude that in the large Nf limit BLM 
is exactly equivalent to z[0/n] PA's. 

In view of this result it is worthwhile to repeat the characteristics of the z[0/n] 
functions, found in Section 4: 



a) As other non-diagonal PA's, a z[0/n] PA does depend on the RS. The leading 
term in this RS dependence is proportional to /3q _1 t n_1 , significantly less than 
pn+i j.n+i - m p ar tial-sum. 

b) As opposed to other non-diagonal PA's, the z[0/n] PA is strictly free from any 
poles in its RS dependence. 



We wish to emphasize that the BLM prescription beyond the leading order (i.e. 
beyond tBLM = ti) cannot be regarded as a choice of RS in the strict sense, since 
terms that depend on z in ( f49"D break the additivity of scale-shifts in the transfor- 
mation (|A.6|) . Putting it differently: after substituting t B LM in (|A.6|) , w depends on 



"This consensus does not include the method described in where tBLM = t\ at any order. 
In this method the higher order terms in {fiox) as well as effects due to sub-leading running 
of the coupling are resummed by setting the higher orders of the (3 function, i.e. choosing the 
renormalization scheme such that the remaining coefficients reach their conformal-limit value. 
This idea cannot be applied in the large Nf limit, where one neglects the higher order corrections 
to the (3 function altogether. 
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z in a more complex way than implied by the RS transformation. This is why the 
BLM result is not free of RS dependence and is not equivalent to a diagonal PA PI. 
This does not apply to the simplest case of a next-to-leading order series, where 
tBLM = t\. Here, BLM in the large Nf limit becomes simply a a;[0/l] PA, which is 
RS invariant. 

We wish to stress that the equivalence of x[0/n] PA's and BLM is of course true 
only in the large Nf limit. Moreover, x[0/n] PA's are different from all the various 
generalization of BLM that attempt to take into account effects due to non-leading 
running of the coupling-constant. If one wishes, one can regard x[0/n] PA's as 
another generalization of this kind. 

7 Conclusions 

We showed that in the (3q limit diagonal PA's of perturbative series become exactly 
renormalization scale independent. 

This implies that diagonal PA's are correctly resumming contributions from 
higher order diagrams which are responsible for the renormalization of the coupling- 
constant. 

Non-diagonal PA's are not exactly invariant even in the large (3q limit, but still 
reduce the global RS dependence as compared to partial-sums. Among the different 
non-diagonal PA's, the only one that has a completely regular behavior with respect 
to scale variations is the x[0/n] PA. We have shown that in the large Nf limit of 
QCD, the latter is identical to the BLM scale-setting procedure. 

In physical cases, higher order corrections in the j3 function break the RS in- 
dependence of PA's, introducing a small scale and scheme dependence, even for 
diagonal PA's. 

We also showed that PA's, when they are indeed RS invariant to a good approx- 
imation, lead to the same numerical result as the ECH method. 

An important feature of PA's (one that in our view makes them more useful than 
scale and scheme setting methods) is simply that they can be used in any scale and 
scheme. The comparison of results in different scales and schemes can then serve 
two goals: 

a) Estimate the reliability of the PA method is in each particular case, by con- 
sidering the scale and scheme dependence of the partial sum as a reference, as 
was done for the Bjorken sum rule in Ref. [Q. 

**The argument we used in subsection 6.2 to show the agreement between ECH and PA's does 
not apply in the case of BLM in the large Nf limit, even though the series reduces to a single term. 
The reason is precisely the fact that BLM is not strictly a choice of RS. 
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b) Use the residual scale and scheme dependence as a lower bound for the theo- 
retical error. 

Finally, we feel that the "surprising success" of PA's in QCD, and generally in 
QFT, is now based on a much more firm basis. 
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Appendix - The large Nf limit 

In this Appendix we briefly present some of the basic formulae that are used 
in the large Nf limit calculations. This concerns the discussion in section 6 and 
especially subsection 6.3 where we compare PA's and BLM in this limit. 

Starting with (Q) and using the fact that the leading term in is proportional 
to Nf 1 , we obtain: 

S(x) = x (l + r\x + r 2 x 2 + r 3 x 3 + • • • + r n x n ^j (A.l) 
= x (l + c x N f x + c 2 N 2 f x 2 + c 3 N 3 f x 3 + ■■■ + c n N]x n ) 

We define z = xN f and S N ^°°(z) = S(x)N f and thus: 

S N r*°°(z) = z (l + Cl z + c 2 z 2 + c 3 z 3 + ■■■ + c n z n ) (A.2) 



In the same manner we change the notations for the /3 function (|^): 

f*£_ Po_ 2 , A_ 3 , & 4 , ... /a ox 



Remembering that /3q ~ A 7 / while at higher orders $ ~ Nf (cf. (f|) and (H)), 
we conclude that higher order corrections to the /3 function are negligible, being 
sub- leading in Nf, and ( |A.3|) can be written as: 

§ = /^°V (A.4) 



where /3 / ^°° = §■ Therefore Eq. translates into: 

z = w + p^°°tw 2 + (p^°°) 2 t 2 w 3 + (p^°°yt 3 w 4 + ■■■ (A.5) 

where w is the coupling-constant that is defined at the new renormalization point \x 
(see Section 2). Finally, Eq. flIU|) translates into: 

(A.6) 



1 _ (3^°°tw 
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RS Transformation for 3 flavours vs. leading fi approximation 
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Figure 1: y(t) as a function of the scale-shift for QCD with 3 flavours, as calculated in 
three different ways: (a) all-order resummation of the leading /3 terms (continuous 
line), (b) first four terms in the leading f3 approximation (dashed line), and (c) the 
first four terms in the actual QCD scale-shift transformation (dotted line). 
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RS Trans, for N f =3 — leading (3 approximation vs. QCD 
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Figure 2: The relative deviation of the leading (3q approximations for y(t) from 
the QCD transformation, given by fl4"6|) . The continuous line represents the rela- 
tive deviation of the all-order resummation of the leading (3q terms (fHj), while the 
dashed line stands for the relative deviation of the first four terms in the leading (3 
transformation (fl5|). 
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